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Abstract: The Q-deformation is a harmonic trap, penning certain excitations near 
the origin in a manner consistent with supersymmetry. Here we explore the dynamics of 
BPS monopoles and vortices in such a trap. We pay particular attention to monopoles 
in the Higgs phase, when they are confined to a vortex string. Unusually for BPS 
solitons, the mass of these confined monopoles is quadratic in the topological charges. 
We compute an index theorem to determine the number of collective coordinates of 
confined monopoles. Despite being restricted to move on a line, we find that they have 
a rich dynamics. As the strength of the trap increases, the number of collective co¬ 
ordinates can change, sometimes with constituent monopoles disappearing, sometimes 
with new ones emerging. 
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1. Introduction 

The O-background is a deformation of supersymmetric gauge theories which breaks 
translational invariance. It has proven to be both a powerful tool for computation 
[1, 2], and a useful device to highlight connections between different theories, most 
notably four-dimensional gauge theories and two dimensional integrable systems [3, 4]. 

In this paper we take a more prosaic view of the O-background. We view it simply as 
a harmonic trap, analogous to those which arise in condensed matter physics. Its role 
is to restrict certain excitations to lie close to the origin. The excitations that we will 
be interested in are solitons. In supersymmetric gauge theories, BPS solitons typically 
have a number of nice properties, both physical and mathematical. The O-background 
provides a harmonic trap which is consistent with supersymmetry and, correspondingly, 
preserves many of these nice properties. 
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There are at least two motivations to study solitons in the ^-background. The first 
is purely classical. A harmonic trap squeezes solitons together. Yet this is often where 
solitons are at their most interesting. They no longer appear as a point-like objects 
and their extended, non-linear nature becomes apparent. They lose their individuality, 
merging into each other to form something new, often with interesting structures and 
collective excitations. 

The second motivation is more quantum in origin. Solitons provide a semi-classical 
springboard to study some of the interesting dualities that are induced by the fi- 
background. In particular, we have in mind the 4d/2d duality described in [5, 6, 7], 
relating the Seiberg-Witten curve to twisted superpotentials of 2d sigma-models. This 
is an extension of an earlier duality [8, 9] which found an explanation in the dynamics 
of vortex strings [10, 11]. Here we study the vortex strings relevant for the extended 
duality. More recently, there have been studies of 3d gauge theories [12] and 5d gauge 
theories [13] in the presence of the fl-background and we will describe the vortices and 
monopoles relevant for these theories. 

What We Do 

We study solitons in J\f = 2 four-dimensional theories with an fl-deformation in a single 
plane. We are not the first to study solitons in this background. A number of properties 
of monopoles were explored in [14, 15] and a range of other solitons were described in 
[16], 

In this paper, we will describe three different types of solitons: monopoles, vortices 
and confined monopoles. Our interest is primarily in the latter. These are monopoles 
in the Higgs phase, where they appear as beads threaded on vortex strings, yet remain 
BPS [17]. Our main results are the formula (2.3) for the mass of the confined monopole 
in a trap, and the index theorem (3.8) for these objects. 

We find that the presence of the harmonic trap endows these confined monopoles 
with a rich dynamics. A generic, higher-charge monopole can split into constituent 
monopoles, each free to move up and down along the string. However, the mass of each 
of these constituents has an extra term which, unusually for BPS solitons, is quadratic 
in the magnetic flux charges. This can be thought of a binding energy between the 
monopole and other flux tubes which also lie in the trap. 

As one increases the strength of the harmonic trap, the number of collective coordi¬ 
nates jumps. Sometimes this reflects the fact that some of the constituent monopoles 
become massless and disappear; sometimes it reflects the fact the new constituent 
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monopoles appear. One of the surprising features is that monopoles with charges that 
one might naively have thought of as anti-BPS can apparently become BPS in the 
presence of the trap. Much of the paper is devoted to telling this story. 

The paper also includes a number of other results. In particular, in two appendices 
we study the dynamics of vortices in the presence of a harmonic trap. The effect of 
the trap is to induce a potential on the vortex moduli space, so that the ground state 
of vortices is an incompressible disc lying at the origin of the plane. We show that, for 
U{1) vortices, the collective excitations of this disc have a description as a field theory 
living on the edge of the disc. 

2. Solitons in a Harmonic Trap 

The theory of interest consists of a U(N ) gauge field A^, coupled to a real adjoint scalar 
c f> and Nf fundamental scalars q t . 

We impose on these fields an external, harmonic trap whose strength is parameterised 
by uj. This is the Id-deformation. The trap breaks translational symmetry, penning 
certain excitations close to the origin in the (a: 1 , a; 2 ) plane; however they remain free to 
move in the x 3 direction. This is the form of the id-deformation discussed in [3]. 

The Lagrangian for this supersymmetric trap was first derived in [18]. (The reference 
[16] contains a useful review of this work.) It is given by 

C = F ' J, ' FIJU ~ ^ Tr W ~ ^ x2El ~ X ' E2))2 - ^ Tr ^ “ “x 1 ^) 2 

N f 

- ^ Tr (^20 - ujx 2 B 3 ) 2 - ^Tr (r> 3 0 + ix(x l Bi + x 2 B 2 )) 2 + 

i— 1 

e 2 Nf 

+ yTr (^ - v 2 ) 2 + |(0 - - iu>(x l V 2 qi + a: 2 £> 1 %)| 2 (2.1) 

i i =1 

Here Vcp = d<fi — i[A, 0] while Vq = dq — iAq. The electric field is Ei = F 0l and the 
magnetic field is Bi = \etjkFjk . In this paper, we focus on the theory with Nf = N ; 
this is the minimal theory admitting vortices. 

The Lagrangian admits a completion to a theory with AT = 2 supersymmetry. In ad¬ 
dition to the fermions, the theory with AT = 2 supersymmetry has a further real adjoint 
scalar in the vector multiplet and Nf anti-fundamental scalars in the hypermultiplets. 
There are non-trivial terms in the Id-deformation involving all these fields [18]. (In 
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particular, there is an extra contribution to the D-term involving the imaginary part of 
the adjoint scalar.) However, it turns out that these extra bosonic fields vanish on our 
soliton backgrounds and, to avoid cluttering equations, we have chosen to omit them 
from the start. 

The Lagrangian (2.1) admits a number of different soliton solutions depending on the 
values of v 2 and the real-valued masses rrq. We now review these different solutions. 

2.1 Monopoles 

When v 2 = 0, we can turn off the fundamental scalars, q t = 0. There is no potential 
for the adjoint scalar 0 and we are free to specify a vacuum expectation value that lies 
in the Cartan sub-algebra. This breaks the gauge group U(N ) —>• U(1) N . 

With these boundary conditions, the theory admits BPS magnetic monopoles. These 
obey the deformed Bogomolnyi equations [14, 15], 

T>](f) — uix 1 B 3 = ±£>i 

V 2 (j) - lox 2 B 3 = ±B 2 (2.2) 

V 3 (j) + u(x 1 B 1 + x 2 B 2 ) = ±B 3 

where the ± signs are for monopoles and anti-monopoles respectively. Solutions to 
these equations describe monopoles carrying magnetic charge in U( l)^ -1 C SU(N), 
trapped at the origin of the (a; 1 , a; 2 ) plane. These monopoles have mass given by the 
usual expression 

M = —J d 3 x d a Tr (f>B a (2.3) 

with a = 1, 2, 3. A study of the simplest, charge one monopole in SU (2) was performed 
in [15]. It was found that the gauge field profile is unchanged by the presence of the 
trap, while the scalar profile of the solution is deformed. In particular, in the presence 
of a monopole, the expectation value of the scalar field (j) differs in different asymptotic 
parts of space. If we pick an expectation value for <f> e su( 2), 

(f> —* a a 3 as a: 1 , x 2 —* oo with x 3 = 0 

then the monopole solution of [15] has 

4> —> (a ± u) cr 3 as x 3 —> Too (2.4) 

This doesn’t affect the mass given in (2.3), which is independent of uj. 
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Little appears to be known about the solutions to (2.2) for magnetic charge n > 2 
or, indeed, for higher-rank gauge groups. In the absence of a trap, it is known than 
there exists a unique axially symmetric 577(2) monopole for each charge n [19]. For 
n > 2, the profiles look like C=E7>. It is not clear if the BPS equations (2.2) have a 
corresponding single solution for charge n > 2 or whether they admit more general 
solutions in which n charge one monopolcs are free to roam in the x 3 direction. 

2.2 Vortices 

When v 2 ^ 0, the fundamental scalar fields q t pick up an expectation value. When 
Nf = N, this is given by 

(Qi) = vS? (2.5) 

The U (. N ) gauge symmetry is now completely broken. In this phase, the theory admits 
BPS non-Abelian vortex strings, first introduced in [20, 21]. 

In general, the adjoint scalar field 0 also picks up an expectation value 

(0) = diag(mi,..., m N ) (2.6) 

The dynamics of vortices depends on these mass parameters m,. We will start by 
describing the vortices when rri t = 0. In Section 2.3 we describe vortices with m, ^ rrij. 

When rrii = 0, the theory has an SU ( N ) global symmetry, in addition to the U ( N ) 
gauge symmetry. The expectation value (2.5) induces the symmetry breaking 

f/(iV) g auge X S17(iV) global —> 517(iV) diag 

This surviving SU (V) diag symmetry acts on the vortices, endowing them with internal 
non-Abelian degrees of freedom. 

In the presence of the trap, the vortex strings are only BPS when aligned along the 
a: 3 direction. The profile of the fields in the (a; 1 , a; 2 ) plane obey the vortex equations 

V z qi = 0 , B 3 = e 2 (q,ql - v 2 ) (2.7) 

where z = x l + ix 2 , together with the supplementary equations for 0, 

V z (/>= l ^zB 3 , d)q t = uz V z q t (2.8) 

Solutions to these equations have tension 

T = —v 2 f d 2 x Tr B 3 = 2irv 2 k k G Z + (2.9) 
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The situation with vortices is rather more straightforward than the situation with 
monopoles. The vortex equations (2.7) remain unchanged which means, in turn, that 
the profiles for A z and q remain the same, ffowever, we now must also solve the equa¬ 
tions (2.8) for <fi in the background of this vortex. These have a simple interpretation: 
they are, respectively, the requirement that A z and q are rotationally invariant, up to 
a gauge transformation. 

The upshot is that, in the presence of the trap, only axially symmetric vortex config¬ 
urations survive. We flesh out this statement in Appendix A, where we also show that 
the dynamics of vortices can be described as motion on the original moduli space in 
the presence of a potential. This potential is generated by the Killing vector associated 
to planar rotations. 

For 7/(1) vortices, there is a unique axially symmetric configuration for each mag¬ 
netic flux k. The zeros of the Higgs field q all sit at the origin. The magnetic flux is 
roughly constant with B « —e 2 v 2 over a disc of radius R ~ \fk/ev , before dropping 
exponentially quickly to zero. Solving (2.8) shows that 0 rises approximately quadrat- 
ically within this disc, <fi ~ uB\z\ 2 /2, before it too drops exponentially quickly to zero 
at a radius \z\ ~ R. 

For vortices in U(N), the situation is more complicated. Each vortex has an internal 
orientation which, when well separated from other vortices, is parameterised by CP jV_1 . 
However, in the presence of the trap the vortices are pushed on top of each other and 
these internal orientations mix in a complicated manner. (We will provide a description 
of this in Section 3.) For k — 2 vortices in 7/(2), the moduli space of axially symmetric 
vortices was determined in [22, 23, 24] to be topologically, but not metrically, CP 2 /Z 2 . 
Surprisingly, this sub-manifold is singular even though the full moduli space of vortices 
is smooth. 

2.3 Confined Monopoles 

When v 2 ^ 0 and m, ^ m.j the theory admits BPS monopoles, but these are now 
confined: the magnetic flux leaves the monopoles along the x 3 direction as a vortex 
string. The resulting solitons are vortex/monopole composites. The possibility of such 
monopoles confined on Z 2 strings was pointed out in [25]. Here we will be interested 
in BPS monopoles confined on non-Abelian strings. These were introduced in [17] and 
further explored in [10, 11, 26, 27, 28, 29, 31, 32, 33]. 

To describe these confined monopoles, it is perhaps simplest to first look at what 
becomes of the vortex strings. With the masses m, turned on, the symmetry breaking 
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pattern is now 


U(N) gmge x SU(N) glob „ ^4 U(N) gme , x {/(I)"!, 1 ,, A U(l)”£ 


where the first of these symmetry breakings is explicit, the second spontaneous. There 
is no longer a surviving non-Abelian symmetry and, correspondingly, the internal orien¬ 
tation modes of the vortex are lifted. Instead, the vortices sit in the Cartan subalgebra 
U(1) N C U(N). The magnetic flux A; of a vortex is distributed among these N Cartan 
elements, and the most general solution takes the form 

B 3 = diag(Bi kl \Bi k *\...,Bl kN) ) and q? = 8?q^ (2.10) 

where 'Yh a k a = k and (B^ k \q^) is the solution for an Abelian vortex with magnetic 
flux k. As we have seen, in the presence of the trap there is a unique Abelian vortex 
solution with a given magnetic flux. This means that the non-Abelian vortex ground 
states are labelled by the ordered set {k a }. 


The number of ground states for a non-Abelian vortex is the number of ways of 
distributing the k total fluxes among the N Cartan elements. In other words, 


Number of ground states 


{k + N-l)\ 
k\(N — 1)! 


( 2 . 11 ) 


The confined monopoles are objects which interpolate from one of these ground states 
to another. They act as sources and sinks for magnetic flux. Because the magnetic 
charge lies in U( l)^ -1 C SU(N ), the monopoles do not change the total flux k. They 
merely redistribute it among the Cartan elements. 


In the absence of the harmonic trap, the Bogomolnyi equations for confined monopoles 
were derived in [17]. In the presence of a trap, it is straightforward to derive modified 
Bogomolnyi equations giving solutions to the equations of motion of (2.1). They are 


T> i<f) — l ox l B - 3 = ±B i 
V-2 <f> — lox 2 B 3 = ±B 2 

V 3 <j) + u(x 1 B 1 + x 2 B 2 ) = ±{B 3 - e 2 { qi q\ - v 2 )) (2.12) 

V z qi = 0 

V 3qi = ±((0 - mi)qi - uzV *&) 


These equations describe BPS vortices (rather than anti-vortices) threaded by either 
BPS or anti-BPS monopoles, depending on the ± sign on the right-hand side. Solutions 
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to these equations have string tension (2.9), together with a finite contribution from the 
mass of the monopole. We will show that this mass is given in terms of the magnetic 
fluxes {k a } at x 3 = ±oo by 


M„ 



m„k n — —k 2 


aj 3 =+oo 
x ^=—oo 


(2.13) 


The first of these terms coincides with the monopole mass (2.3) when cj) has expectation 
value (2.6), albeit with the flux now collected at x 3 = ±oo as befits a vortex string, 
rather than radially for a monopole in the deconfined phase. The second term, quadratic 
in fluxes, is novel. As we will explain shortly, it describes the binding energy between 
the monopole and vortex string. 


Computing the Mass of the Confined Monopole 

We now derive the formula (2.13). We look for static solutions with d 0 = A 0 = 0. We 
take the energy functional associated to (2.1) and complete the square thus: 


£ = 


f ^ X 2e2 Tl ' _ ux'B-s + eBi) 2 + ^Tr (V 2 <j) - ujx 2 B 3 + eB 2 )‘ 


1 

2e^ 

N 


Tr 


V 3 (j) + u>x 1 Bi + ujx 2 B 2 + e(B 3 - e 2 (^ g^} - v 2 )) 


i =1 


+ ^ < | T> z qi\ 2 + \V 3 q t - e((0 - rn^qi - iuj(x l V 2 qi - x 2 !)^))^ 


—v 2 Tv B 3 — e—Tr B a V a <t> — ev 2 TrV 3 (f> — eujv 2 Tr (x l B 3 + x 2 B 2 


N 


1=1 


eV 3 (qj(cf) - rm)qi ) + ec oqj(x 1 B 1 + x 2 B 2 )qi 


-eiujV 3 q\(x l V 2 qi - x 2 T>iqi) + eiu)(x L V 2 q} - x 2 T>iq})V 3 q, 


/nft ) 


where e = ± determines whether we’re dealing with monopoles or anti-monopoles. The 
Bogomolnyi equations (2.12) can be seen in the total squares of the first three lines, 
while the first term of the fourth line is the tension of the vortex string. The remaining 
terms are the mass of the confined monopole. After integrating by parts, and using 
T> a B a = 0, we can write them as 

Mmono = e f d 3 x - -^<9 Q: Tr (B a (j>) - v 2 TrV 3 (j) - ujv 2 Tr (a^Bi + x 2 B 2 ) 



(2.14) 


+ 



All these terms are total derivatives. As we already mentioned, because the flux is col¬ 
limated in a vortex string, the first term gets contributions from x 3 = ±oo rather than 
radially. The last of these terms vanishes because, asymptotically, the configuration 
looks like an axially symmetric vortex string, with 0 obeying (2.8). We’re left having 
to evaluate 


^mono ^ ^2 I d X 


]2 ~ Tr 4>(B 3 + e 2 v 2 ) + coe 2 v 2 Ti (—x 1 A 2 + x 2 Ai) 
We can massage the last term. We write 


x 3 =+oo 

x 3 =—oo 


UJ 


d 2 x Tr (— x 1 A 2 + x 2 Ai) = u 


d 2 X yTl'5,3 


d 2 x Tr 0 


where we’ve used the expression ( 2 . 8 ) for the profile of 0 in the background of the 
vortex. We’ve also discarded some boundary terms upon integration by parts; these 
are common to a : 3 = ±oo and so do not contribute to the monopole mass, which 
becomes 


A/mono T‘ 


d~x 


2 „. 2 \ 


Tr 4>(B 3 + 2e v 


aw=+oo 


x°=—oo 


(2.15) 


This is pleasingly simple, involving asymptotic integrals at x 3 = ±oo. However, it’s 
not obvious that it depends only on topological information. Indeed, this is even true 
of the first term which, naively, looks the same as the usual monopole mass (2.3). The 
reason for this is that, as we saw in Section 2.2, the field 0 gets a profile in the presence 
of a vortex. This profile is complicated, determined by the equations (2.8), and feeds 
into the expression for the mass. Nonetheless, as we now show, nice things happen. 


The nice things follow from the fact that the asymptotic vortices are axially symmet¬ 
ric, so that both \q\ 2 and 0 are a function only of the radial coordinate r. Consequently, 
it is straightforward to show 


<9 r 0 = uirB-i 


with 0 —> ( 0 ) = diag(mi,... , tttat) as r —> oo, and 0 —* ( 0 ) + cadiag(£;i,..., k^) as 
r —> 0, a result which was also used in [16]. Using all these facts, we have 


J d 2 x Tr 0 I ?3 


2vr 



-2vr (fc*™* + 7 \ k f) 

i 


(2.16) 


We’re left with having to compute the linear term, f d 2 x Tr 0 in (2.15). This follows 
straightforwardly from the vortex equations (2.7) and (2.8), which imply 
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Tr 0(83 + eV) = e 2 g?0® 

i 

= e 2 ^ ////<'/•<'// + uqjzV.qi 

i 

= "d (-83 } + e2 ^ 2 ) + ue 2 V- z (z(q\qi - v 2 )) - cue 2 (gjg; - v 2 ) 

i 

= ^ m, ^ 83 ^ + e 2 u 2 j — a; 8 g l ' ) 

i 


where the magnetic held has been decomposed into its Abelian components (2.10). The 
integrals can now be expressed in terms of the Abelian fluxes {k a }. Dropping terms 
common to the two boundaries at x 3 = ± 00 , we have 


/ 


d 2 x 


Tr 0(83 + e 2 v 2 


£ 3 =+00 
x 3 =—oc 



Comparing to (2.16), we learn that 


e 2 v 2 / d 2 x 


Tr, 


x a =+oo --, 

3 =+ 27 r 5 Z 

x°=—oo 


\-k 2 
. 2 1 


£ 3 =+00 

x 3 =—oo 


Finally, we can put this all together to give 


27r ^ 
h i mono = =F 2 _^ 


m,ki - z ir 


(2.17) 


x = —00 


which is the promised result (2.13). 


Interpreting the Monopole Mass 

The term linear in k is the familiar monopolc mass term. However, it is unusual to find 
BPS solitons whose mass is quadratic in the charges. Here we provide an interpretation 
of the quadratic term in (2.17). We will argue that it can be thought of as a binding 
energy between the confined monopolc and vortex strings. 

To see this, it will suffice to work with U (2) gauge theory. We choose the masses to 
be mi = — m and m 2 = +m, so that the asymptotic expectation value is ( 0 ) = — ma 3 . 
Let’s first look at the charge k monopole, interpolating between the vortex strings with 
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flux (ki, k-i) = (0, k ) at x 3 = — oo and flux (fc, 0) at x 3 = +oo. The confined monopole 
looks like: 


B = 



o m ) 


The mass of the monopole is given by 


Mmono 


Anmk 


(2.18) 


This is the usual mass of a monopole in an SU{2) gauge theory. There is no contribution 
from the harmonic trap. 


Let’s now look at the same monopole of charge k , this time interpolating between 
(j)i,k + p 2 ) at x 3 = —00 and (k + pi,p 2 ) at x 3 = + 00 . We can think of this as a 
combination of the previous configuration, superposed with a vortex with flux (pi,p 2 ). 

.) ( 

I 


(I 


The mass of the monopole is now 

4:7imk 2nu)k 

^mono I - ~^2 vPl P 2/ 

The second term is due to the presence of the (pi,p 2 ) vortex string and, in this sense, 
can be thought of as the binding energy between the monopole and this string. Note 
that this binding energy can be positive or negative. Note also that oj can take either 
sign. 
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There are further questions that we would like to ask. Let’s return to the charge 
k monopole with mass (2.18). Can this monopole decompose into charge one con¬ 
stituents? The first few monopoles in the string would look like this: 



If this is possible, each charge one monopole now experiences a binding energy with the 
underlying vortex string. The mass of the q th monopole in the sequence (q — 1 ,k) 
is given by 


M q 


4:71m 27 TUJ 


{k — 2q + 1) 


(2.19) 


Although the mass of each individual constituent gets a contribution from the binding 
energy, the sum of all these masses is again (2.18), independent of oo. 


There’s something unsettling about the mass formula (2.19). For the sake of the 
argument, let’s choose cu > 0. The first monopoles in the sequence (2.19) then have 
negative binding energy; the later ones positive. More worryingly, if we take this 
formula at face value then it looks as if the mass of the first monopoles will become 
negative for sufficiently large k. This is clearly unacceptable. 


In the next section we will understand what’s going on with this necklace of monopoles. 
We will show that a charge k monopole can, in fact, decompose into its constituent 
parts as described above, but this can only happen if the masses of those constituents 
are positive. Moreover, the story for higher rank gauge groups will turn out to be 
somewhat more involved. 


3. Monopoles as Kinks on the Worldsheet 

In this section, we explore the confined monopoles in more detail. However, instead of 
solving the Bogomolnyi equations (2.12), we will instead look at the monopoles from 
the perspective of the vortex string worldsheet, where they appear as BPS kinks. This 
was first shown in [17], and further explored in [10, 11, 30]. However, attention has 
previously focussed on monopoles threaded on a single k — 1 vortex string. These are 
restricted to a single monopole in 1/(2) gauge theories or the so-called g = (1,1,..., 1) 
monopoles in U(N ) theories. 
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Here we will be interested in monopoles threaded on more than one vortex string. In 
part, the problem with studying this in the past was to decouple the effect of the string 
positions from the dynamics of the monopoles along in the string. Here, the harmonic 
trap helps us since, as we have seen, the strings are obliged to sit at the origin. 

3.1 A Matrix Model for Vortex Strings 

We could study the vortex string dynamics from the geometric perspective of the moduli 
space, as described in Appendix A. However, a more versatile description is offered by 
the vortex matrix model introduced in [20]. 

The matrix model provides a description of the vortex dynamics in terms of an 
Af = (2, 2) gauged linear sigma model. It is a U(k) gauge theory, with a single complex 
adjoint scalar Z and N fundamental scalars ipi, i = 1,..., N. We also include an 
auxiliary adjoint scalar, a. Neglecting fermions, the low-energy vortex dynamics is 
then described by 

N n 2 

Tr \V Q Z\ 2 + Y, I V a <pi\ 2 - yTr ([Z, Z t] + ^ - r) 2 

i=l i 

l N 

+ X 2 Tl 0 D a<?) 2 - Tr |[<7, Z\ + UjZ I 2 - Y <Pi( a ~ l 3 - 1 ) 

& i=l 

Here a = 0, 3 runs over the vortex string worldsheet directions. Note that although a 
is generally a complex adjoint scalar, in the present case where we have restricted to 
real masses rn t , we can take a to be real. 

We should think of this action in the g 2 —y oo limit. (There is also a kinetic term 
for the U(k) gauge field which we have omitted; it is suppressed in this limit.) The 
“D-term” potential on the first line is then imposed as a k x k matrix constraint 

N 

\ Z ' Zt ] + VMi = rlk ( 3 - 2 ) 

i= 1 

If we further quotient by the U(k) gauge action, Z —> UZW and ip —>• Uip, this defines 
the vortex moduli space A4 It has dimension 

dim(Al N,k) — 2 Nk 

Note that this quotient construction defines a metric on M.N,k- This does not coincide 
with the metric on the vortex moduli space defined in Appendix A. Nonetheless, this 
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construction will suffice for our purposes. In particular, the masses of BPS states 
coincide with those in the vortex theory if we fix the Kahler class of M. N.k to be [20] 


r = 


2vr 


The auxiliary scalar a can be ignored if to — = 0. In this case, after the quotient 

above, the action (3.1) becomes a non-linear sigma model on M.N,k- However, if either 
u 0 or rrii ^ 0 then a is non-vanishing and the last two terms in (3.1) can be thought 
of as inducing a potential on M.N,k- 


U(l) Vortices 

Let us look at vortices in the U( 1) theory in more detail. Here there is a single p held 
in the action (3.1). We may set the mass of p to zero through a constant shift of a. 
The ground state of k vortices is then given by the constraint (3.2), supplemented by 


(7 , Z] + uZ = 0 


(3.3) 


There is a unique solution to these two equations, found in [48]. This is 


Zk 


(° 1 

0 y/2 


y/r 


\ 


and p 


/ 0 \ 

I 0 


y/r 


(3.4) 


with a = a k = —uj diag(h — 1, k — 2,..,, 2,1, 0). 
rotationally invariant k-v ortex solution. 


This corresponds to the unique, 


U(N) Vortices 

For U(N) vortices, the situation depends on the masses m,. For rrp = 0 (or, indeed, 
m, all equal) there remains a moduli space of solutions to (3.2) and (3.3). 

Here, we are instead interested in the case of distinct masses rn t ^ rrij. There are 
now isolated solutions to (3.2) and (3.3) in which the total flux k is distributed among 
the N Cartan elements, mimicking our discussion in Section 2.3. The general solution 
is determined by an ordered set {ki} with JV /q = k. The matrix Z is block diagonal, 
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while o is diagonal, although conveniently written in block diagonal form 


( Z, 


Z = 


k\ 


\ 


J k2 


Zk » 


/ 


and a = 


a kl + m i 


\ 


Vk 2 + m 2 


(3.5) 


+ m N J 


with Zfc given by (3.4) and ay by the expression below (3.4). Similarly, each yy has only 
one non-zero element, y/rkl. The number of solutions of this form is given by (2.11). 


3.2 Kinks on the Worldsheet 

When u ^ 0 and rn l ^ rrij , the worldsheet theory has a large number of isolated ground 
states. There exist BPS kinks which interpolate between these some of these ground 
states. These obey the Bogomolnyi equations 

V-itpi = (a - ny)^ 

V 3 Z = [a, Z] + uZ (3.6) 

V 3 a = g 2 ([Z, Zt] + ipiipj - r) 

BPS domain wall equations with a similar structure have been studied previously, both 
in Abelian [34, 35] and in non-Abelian [36] gauge theories. The novelty here is the 
presence of the adjoint scalar Z. As we will see, this changes the dynamics of these 
kinks. 


Writing the energy as a sum of squares, it is simple to show that the mass of any 
kink obeying these equations is given by 


Afkink — — 


y] nii(p\(pi — cuTr Z^Z 


-i # 3 =H-oo 


x ,J =—oo 


We can evaluate this on the ground states (3.5). In the vacuum {/q}, we have yd</y = 
rk tl and 

N 

TVZ'Z = 5 - !) = k ‘ - *) 


i =1 


The mass of the kink is therefore given by 

N 


Mkink = -r ^2 


1=1 L 


UJ , 


rrijkj — —ki 


x 6 =-\-oo 


j x°=—oo 


(3.7) 


Identifying r = 27r/e 2 , this agrees with the mass of the monopole (2.3). 
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3.3 Counting Zero Modes for Confined Monopoles 

We now return to the question that we asked in Section 2.3: can a monopole of given 
charge decompose into constituent monopoles, spread like a necklace of beads along 
a vortex string? To answer this question, we attempt to count the number of zero 
modes of the kinks. We do this by computing the index for the linearised domain wall 
equations (3.6). This calculation is presented in Appendix C. The result is: 


X 


1 

2 


ki-l 

EE sign(mj — rrii + (q — kj)u ) 

. i,j <j=0 


# 3 =+oo 


#3 —_qo 


(3.8) 


The sum is over all i, j (i.e. each pair is counted twice). It is understood that there is 
no contribution from any sector where k t — 0. 


We will now illustrate the physics underling this index theorem by giving a few simple 
examples. But first we give some important warnings. The index counts the number of 
zero modes of the operator A arising from linearised Bogomolnyi equations, minus the 
number of zero modes of the adjoint operator Ah (Both of these operators are defined 
in Appendix C.) Usually, when computing the index for BPS solitons, one can show 
that the adjoint operator has no zero modes. This means that the index counts what 
we care about: the number of collective coordinates of the soliton. This is the case, 
for example, in [37, 38, 39, 40] and [29]. It is not, however, the case here. There are 
regimes of the parameters m t and u for which we can show that the adjoint operator 
does have zero modes. This means that the index X can only give a lower bound on 
the number of collective coordinates. 


On top of this, although we can count the number of zero modes of solutions, we have 
no proof that solutions actually exist, at least beyond the simplest cases with Z = 0 
[41]. These issues mean that any attempt to extract physics from (3.8) alone will be far 
from rigorous. Nonetheless, if we assume that solutions to (3.6) exist whenever X > 0, 
then a consistent story emerges from a combination of the index and the mass formula 
(3.7). We now tell this story. 


Confined Monopoles in U( 2) 

We start by revisiting the problem that we introduced in Section 2.3: is it possible for 
a charge k SU{2) monopole to decompose into its constituent charge one components? 
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The first few monopoles of the resulting configuration would look like this: 


k- 1 


k- 2 


k — 3 


III 


Thankfully, for U (2) monopoles, or domain walls with two flavours, the adjoint operator 
At has no zero modes, so we take the index X at face value as counting the number of 
collective coordinates of our confined monopole. 

As in Section 2.3, we set m* = (— m, +m) and take cu > 0. The index theorem (3.8) 
becomes 

1 k 1 

X = - [sign(2m + qoj) — sign(—2m + qco)\ = -[sign(—2m + qu)\ 

q=0 q=0 

We see that the number of zero modes depends on the relative values of m and c o. In 
what follows, we fix the masses and see how the index varies as we change oo. Suppose 
that we start with small oo, so the trap is weak. Then we have 

(k — l)u; < 2 m =>- X = k 

This counts complex zero modes. It strongly suggests that the charge k monopole 
can indeed split into its constituents. Each of these constituents has a position, and a 
phase. The phase associated to domain walls was first discussed in [42] and is entirely 
analogous to the phases associated to monopoles; in both cases, excitations of this 
phase coordinate turn the monopoles into dyons. 

However, as we increase the strength of the trap, interesting things happen. We 
see from (2.19) that at the threshold (k — l)cu = 2m the mass of the first monopole 
vanishes: Mi = 0. Moreover, the index theorem (3.8) tell us that, at this point, the 
number of zero modes jumps by one: 

(k — 2 )oj < 2m < (k — l)o; X = k — 1 

What’s happening here is the following: as (k — 1)cj —> 2m, the left-hand monopole 
becomes massless and swells to infinite size, joining with the second monopole on the 
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string. The first few monopoles are now 


k- 2 


k -3 




The first monopole on the string now has charge 2, and the others charge 1. From 
(2.13), the mass of this charge 2 monopole is 

AnT 

M' = — (2 m — l oik — 2)) 
e- 

The masses of the subsequent monopoles are still given by (2.19). As we lower the mass 
still further, we again hit a threshold where this charge 2 monopole becomes massless 
and the index (3.8) jumps once again 

(k — 3 )oj < 2m < (k — 2)c a =>• 1 = k — 2 

This pattern now continues. The necklace consists of a charge 3 monopole followed by 
a string of charge 1 monopoles. As the mass decreases, it is always the left-most, higher 
charge monopole which becomes massless. In general, if the mass lies in the range 

(k — p — l)(u < 2m < (k — p)uj =>- 1 = k — p 

there are k — p constituent objects, which arrange themselves on the string as a single 
charge p + 1 object, followed by followed by k — p — 1 charge one monopoles. By the 
time we get to oj > 2m, there is just a single zero mode left. 

We learn that a monopole can split into constituent parts only if the mass of these 
constituents is positive. As we increase the strength of the harmonic trap, the quadratic 
binding energy causes more and more monopoles to become massless and there are fewer 
physical constituents. It seems plausible that the jump of the index at integer values of 
2 m/uj is related to the shift of the asymptotic value of <j) (2.4) found in [15]. It would 
be interesting to understand this connection better. 

Confined Monopoles in U( 3) 

We now turn to monopoles in 17(3). Here, for the first time, we run into the issue of 
zero modes for the adjoint operator A^, which means that the index (3.8) undercounts 
the number of zero modes. For this reason, we go slowly and present three examples of 
successive complexity. In all of these, we again order the masses such that rn l+ \ > m l . 
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Example 1: 


We start by considering the example of flux {ki} = (0,1,1) at x 3 —oo interpolating 
to {kj,} = (1, 0,1) at a; 3 = +oo. We write this as 

( 0 , 1 , 1 ) —► ( 1 , 0 , 1 ) 

This is simply a f/(2) monopolc embedded in the upper-left corner of 1/(3), now ac¬ 
companied by a single flux tube in the lower-right corner. One might think that these 
two objects don’t talk to each other and, indeed, for most values of u the index is 
X = 1. But, for u in a particular range, the index tells us that something interesting 
can happen. We have 

{ 1 u < m 3 — m 2 

1 u > m 3 — mi 

2 otherwise 

In the window (m 3 — m2) < u < (m3 — mi), it appears that this single monopolc can 
decompose into two constituents. We believe that this is indeed the correct interpre¬ 
tation. A study of the asymptotics of fluctuations around (3.6) suggest that, for this 
range of ca, the scalar Z is not just a bystander in the kink equations; its off-diagonal 
components can turn on. This corresponds to a the charge one monopole decomposing 
into the following two constituents: 

(0,1,1) —► (0,0,2) —► (1,0,1) 

At first glance, this looks unlikely. The first monopolc has charge (0, —1,1), the second 
charge (1, 0, —1). But naively (where “naive” means “in the absence of the trap”), the 
second monopole would be BPS and the first anti-BPS. It is very unusual for these 
two monopoles to sit happily next to each other because they carry opposite magnetic 
charge in the final component. (Something similar, although not entirely the same, 
happens with calorons [43, 44].) 

To see why this is acceptable in the harmonic trap, let us look at the mass formula 
(3.7). This is the mass formula for a BPS monopolc. Usually, anti-monopoles would 
have negative mass when evaluated on this formula, reflecting the fact that we should 
use the anti-BPS mass formula which differs by a minus sign. However, we find that 
the mass of the monopolc interpolating from (0,1,1) —> (0, 0, 2) is indeed positive when 
u > (m 3 - m 2 ). 
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The emergence of this new monopole of positive mass explains the jump in the index. 
But why does it subsequently decrease back to I = 1 when uo > (m 3 -mi)? This can be 
seen by looking at the second monopole, interpolating from (0,0,2) —» (1,0,1). From 
(3.7), we learn that this monopole has negative mass when to > (m3 — mi). In this 
regime, we are back down to a single bead on the vortex string. 

Example 2: 


Let us now consider a very similar example, 

( 1 , 0 , 1 ) —► ( 1 , 1 , 0 ) 

This is again a I/(2) monopole accompanied by a I/( 1) flux in a different sector. As 
in our previous example, there is a window of to in which the index exhibits a jump. 
Except, this time it goes the other way: 


{ 1 to < m 2 — mi 

1 to > m 3 - mi 

0 otherwise 

In the window (m2 — mj) < to < (m 3 — mi), the index vanishes. However, we know 
that this cannot be a correct count of the number of zero modes because it is always 
possible to embed the known U (2) solution, with its single complex zero mode, in the 
lower-right corner of the U (3) matrix. 

What’s happening here is that, in this window of to, the adjoint operator gains a 
zero mode. (Again, this conclusion is reached by investigating the asymptotics of the 
appropriate linearised equations.) This is not a physical fluctuation, but is subtracted 
from the index. This means that the index is under counting. For this example, the 
correct number of (complex) collective coordinates of the confined monopole is 1 for all 
values of to. 

Example 3: 


As our final example, we consider a U (3) monopole which, in the absence of a trap, 
would decompose into constituents [45, 39]. This is the monopole 


(0,1,1) 


(1,1,0) 
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The index (3.8) is 


X = 2 + - [sign(m 2 — mi — to) — sign(m 3 — m 2 — a;)] 

For small to, the index is X = 2 suggesting that the monopole can split into two 
constituents. These are naturally identified as the (1, —1, 0) monopole and the (0,1, —1) 
monopole. There are two possible orderings of these constituents on the vortex string, 
resulting in the sequence of flux tubes 

( 0 , 1 , 1 ) —► ( 0 , 2 , 0 ) —► ( 1 , 1 , 0 ) 

(0,1,1) —► (1,0,1) —► (1,1,0) (3.9) 

We conjecture that both orderings are allowed. In other words, there is a moduli space 
of monopoles of these two different types and these two monopoles can pass through 
each other. 

Now we increase to. There is a window where the index either jumps up to 3 or 
down to 1. When the index rises, we would like to understand physically what this 
corresponds to. When the index drops, we would like to understand if this reflects 
a genuine drop in the number of collective coordinates or if it is instead a fake drop 
caused by the emergence of zero modes for the adjoint operator. 

The physics is different depending on the specific values of the masses. There are 
two cases to consider: 


Case A: m 2 — m\ < m 3 — m 2 


In the window, (m 2 — mi) < to < (m 3 — m 2 ), the index decreases to X = 1. Following 
our discussion of U (2) monopoles, we’d like to understand this in terms of one of the 
constituents becoming massless. Indeed, if we look at the first route in the ordering 
(3.9), we see that the monopole (0,2,0) —> (1,1,0) does indeed become massless at 
to = (m 2 — m i) as we enter this window. This means that, inside the window, we 
cannot decompose the monopoles in the first ordering. 

However, none of the monopoles in the second ordering in (3.9) become massless as we 
enter the window. From this, we conjecture that the decrease in the index is actually a 
fake: the number of collective coordinates does not decrease at to = (m 2 —mi). Instead, 
the decrease in the index is telling us that only the second ordering of monopoles is 
allowed. In fact, the two constituent monopoles in the second route remain massive for 
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all values of oj. We conjecture that the constituent structure of monopoles in this case 
is given by 


u j < (m 2 — mi) Two constituents, with either of the orderings (3.9) 

u j > (m 2 — mi) Two constituents, with only the second ordering in (3.9) 


Case B: m 3 — m 2 < m 2 — mi 

In the window, (m 3 — m 2 ) < oj < (m 2 — mi), the index now increases to X = 3. As 
in previous examples, the lower edge of this window reflects the emergence of a new 
monopole with positive mass. There are now three constituents, but there remain only 
two possible orderings because the newly emerged monopole only has positive mass if 
it sits first on the string: 

( 0 , 1 , 1 ) —► ( 0 , 0 , 2 ) —► ( 0 , 2 , 0 ) —► ( 1 , 1 , 0 ) 

(0,1,1) —► (0,0,2) —> (1,0,1) —> (1,1,0) (3.10) 

The first monopole in the sequence has charge (0,-1,+1) and has positive mass for 
oj > ( m 3 - m 2 ). 

We now increase oj further to reach the upper end of the window at oj = (■ m 2 — mi ). If 
we have taken the first route, the final monopole (0,2,0) —y (1,1,0) becomes massless 

at this point and must merge with the middle monopole, leaving us with the sequence 

(0,1,1) —* (0,0,2) —> (1,1,0). However, this is a subsequence of the second route. 
And, from the perspective of the second route, nothing bad happens at the point 
u j = (m 2 — mi). We therefore suspect that this upper end of the window is a fake and 
the number of zero modes does not decrease from 3 back to 2 at oj — (■ m 2 — mi). 

However, this is not to say that the number of zero modes remains 3 for all oj. 
Something interesting does happen at the later point u ; = (m 3 — mi). Here the monopole 
(0, 0, 2) —> (1, 0,1) becomes massless. 

Putting these facts together leaves us with the following conjecture for the decom¬ 
position of monopoles in this system 

oj < (m 3 — m 2 ) Two constituents, with either (3.9) ordering 

( m 3 — m 2 ) < oj < {rn 2 — mi) Three constituents, with either (3.10) ordering 
(m 2 — mi) < oj < (m 3 — mi) Three constituents, with only the second (3.10) ordering 
oj > (m 3 — mi) Two constituents, with only the second (3.9) ordering 

Notice that, for both Case A and Case B, we have found ourselves in the same place 
for large oj, with a restriction on the ordering of the monopoles. 
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The General Story 


The examples described above seem to capture the kinds of behaviour that can arise 
for monopoles in U (. N) gauge groups. In general, we conjecture the following rules 

• If the index increases, it can be traced to a new constituent monopole moving 
from negative to positive mass. 

• If, for all orderings, some constituent monopole becomes massless, then the num¬ 
ber of collective coordinates decreases. 

• If the index decreases and a constituent monopole becomes massless only for some 
orderings, then those orderings are no longer allowed. The number of collective 
coordinates does not decrease. 


4. Summary and Discussion 

We hope that the discussion above has illustrated the rich dynamics of confined monopoles 
in the presence of a harmonic trap. Clearly there is much still to understand. We have 
formulated a conjecture for how monopoles decompose into different constituents with 
different orderings but, because the index is not equal to the number of collective co¬ 
ordinates, we have had to rely on some guesswork and physical intuition. Ideally, this 
guesswork can be placed on firmer footing. 

If our conjectures are correct, then the next question is to get a better handle on the 
moduli space of these confined monopolcs and, in particular, how this moduli space 
changes as we vary oj. 

For a small subset of confined monopoles, the moduli space has been well studied 
[36, 29]. These are the kinks with Z = 0. They only exist when the number of vortices 
is no greater than the rank of the gauge group, k < N, and even then, describe only 
a small number of the possible confined monopoles. It was shown in [30], there is a 
close connection between these Z = 0 kinks and magnetic monopoles in the absence 
of a trap. The kink moduli space was shown to be a middle-dimensional complex 
submanifold of the corresponding monopole moduli space. This submanifold is defined 
as the fixed point of a U( 1) isometry which rotates the monopole configurations in the 
(a; 1 , a; 2 ) plane. 
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One may wonder whether there is a similar connection between free monopoles and 
the more general confined monopoles with Z / 0. If such a connection exists, it must 
be somewhat more complicated. It is known, for example, that the Atiyah-Hit chin 
moduli space, describing charge two SU(2 ) monopoles, has a unique fixed point under 
rotations. Yet we have conjectured above that the corresponding confined monopoles 
have a moduli space in which the constituents can move independently along the string. 
To find a connection in the general case, one should presumably first get a better handle 
on the space of solutions to the monopole equations (2.2) in the presence of trap. 
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A. Appendix: Dynamics on the Vortex Moduli Space 

In this appendix we show that the dynamics of vortices in the presence of the fl- 
deformation is governed by motion on the usual vortex moduli space, but with the 
addition of a potential. This potential is the length-squared of a Killing vector, K, 
associated to action of rotations. 

This result is not surprising and no doubt known to experts. Indeed, the original 
purpose of the O-deformation was to introduce exactly such a potential on the instanton 
moduli space [1]. Nonetheless, to our knowledge an explicit derivation has not appeared 
previously in the literature. Here we closely follows the method of [46] which showed 
how related potentials on other soliton moduli spaces can arise. 

Let us start by reviewing the dynamics of vortices in the absence of a trap. The 
moduli space A4N,k is defined to be the space of solutions to the U(N) vortex equations 
(2.7) with magnetic charge k. Index theorems show that [37, 20] 

dim(.Mjv,fc) = 2 Nk 

The low-energy dynamics of vortices can be described by motion on M.N,k [47]. As we 
now review, the kinetic terms in (2.1) induce a natural metric on M.N,k such that the 
path of followed by vortices as they scatter is a geodesic on M.N,k- 
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We start by introducing coordinates X a , a = 1,. .., 2Nk on M.N,k■ These are the 
collective coordinates of the vortices, meaning that the most general vortex solution 
can be written as A z = A z (x] X ) and qi = q, (x] A"). We define the zero modes to be 

c)y\ do ■ 

ti a A z = + VA-la and 5 a qi = + iVL a qi 

where 12 a (o;; A") is a gauge transformation. By construction, these obey the linearised 
version of the vortex equations (2.7) for any choice of 12 a . To determine 12 a , we need 
to fix a gauge. We ask that the zero modes obey the following equations: 

T> z 8 a qi - i5 a A z qi = 0 and V z 5 a A z = —— 5 a qi qj (A.l) 

The first of these is the linearised version of the first vortex equation in (2.7); it is 
obeyed for any choice of 12. The second equation combines the linearised second vortex 
equation in (2.7) with a gauge fixing condition; it determines 12 a . 

The key idea of the moduli space approximation is that we can accurately describe 
vortex dynamics by simply promoting X a —> X a (t,cr) where a = x 3 parameterises the 
spatial extent of the string. If we further set Aq = —Q a X a and A 3 = —Q a X a ', then the 
kinetic terms for fields are given in terms of the zero modes 

V a qi = 8 a qi d a X a and F ai = 5 a Ai d a X a a = 0, 3 

The gauge fixing condition in (A.l) has been chosen so that these fields satisfy the 
Gauss’ law constraints from the gauge theory. 

Substituting this ansatz into the kinetic terms in (2.1) gives us our expression for 
the vortex dynamics in terms of a non-linear sigma-model on M.N,k, 

S v ortex = I dtda g ab (X ) d a X a d a X b 
where the metric g ab is given by 

r 2 ^ 

gab = I d 2 x — Tr ( S a A z 5 b A z + S a A z S b A z ) + Y]( 8 a ql 8 b qi + 5 b qjd a qi) 

Adding the Harmonic Trap 

Now let’s see how things change when we add the harmonic trap. In general, this 
increases the energy of a vortex solution. When evaluated on the vortex equations 
(2.7), the potential terms in (2.1) simplify and can be written as 

r 2 N 

V = d 2 x — Tr I V z (j) - %-zB 3 \ 2 + ^ - uizV-q ^ 2 (A.2) 

J e 1 i= 1 
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The potential vanishes only if equations (2.8) are obeyed. These are the surviving, 
static vortices in the presence of a trap. For other vortex configurations, we should 
determine 0 by minimising V in the vortex background. In this manner, the potential 
V can be thought of a potential over the moduli space M.N,k- 

It is useful to have an expression for the potential V purely in terms of moduli space 
data. We will prove the following 

C = toA- r VC (A.3) 

where K rot is the Killing vector on M.N,k arising from spatial rotations in the (a; 1 , a; 2 ) 
plane. 

To show this, we first need to examine more closely the zero modes associated to 
spatial rotations. They take the form 

SrotQi = iu(zd z - zd~)qi + ittqi 
<5rot A z = uj(zd z — zd z )A~ + iA z + T) z i} 

where the extra term in 5 ro tA z reflects the fact that the field has spin 1. We’ve chosen 
to normalise the zero modes by u, the strength of the trap. To make progress, we need 
to identify the gauge transformation We will see that it is useful to decompose the 
gauge transformation as 

= — iuj(zA z — zA z ) + 0 

Note that we’ve elected to call this gauge transformation 0. This is not a coincidence. 
With this choice the rotational zero modes can be conveniently written as 

£rot<?i = -izojV z q, + i(j)qi 

5rotA z = -^-B 3 + Vj (A.4) 

It remains to determine 0. As we have described, this is fixed by the second of the 
equations in (A.l) which becomes 



But this is rather nice. It is precisely the equation obeyed by the field 0 that arises from 
minimising (A.2) in the background of a vortex. In other words, 0 = 0. In particular, 
this highlights a claim that we made previously: for rotationally invariant vortices, the 
field 0 can be thought of as a compensating gauge transformation. 
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The overlap of the rotational zero modes gives the length squared of the Killing 
vectors. From (A.4), we see that it can be written as 

, 2 N 

gabK ot KL = / d 2 x —Tr \V z <f> - ^zB 3 | 2 + ^ |0& - ojzV^q,} 2 
d e i= 1 

which proves our claim (A.3). 

In the presence of a harmonic trap, there are also new contributions to the kinetic 
terms, coming from T> 0 (j) terms in (2.1). These are of order O(b 0 2 X 2 )] they can be 
ignored only if we think of both derivatives and the potential being small. If this is the 
case, the moduli space approximation for the dynamics in the presence of the trap is 
given by 

^vortex — I dtda g ah (X) d a X a d a X b - g c jX)K° ot K h rot 

B. Appendix: Dynamics of U( 1) Vortices in a Trap 

In this appendix, we describe the dynamics of 17(1) vortices in a harmonic trap. For 
simplicity, we focus on vortex particles in 7 = 2 + 1 dimensions rather than the vortex 
strings in d — 3 + 1 described elsewhere in this paper. 

Each vortex is a disc of magnetic field of radius ~ 1/ev. As k vortices are brought 
together, they form a larger disc of radius R ~ y/k/ev. This fact has motivated the bag 
model of vortices [49, 50] in which vortices are thought of as an incompressible fluid. 
However, applications of the bag model for BPS vortices have always been limited by 
the fact that vortices could fragment at zero cost of energy. 

The presence of a harmonic trap prohibits this fragmentation. Now we can ask what 
the low-energy dynamics of BPS vortices looks like. We will show that, as expected 
from the bag model, the low-energy dynamics of k vortices can be described by a free 
scalar field propagating around the edge of the disc. 

We work with the matrix model of vortex dynamics described in Section 3. (It would 
be interesting to revisit this calculation using a more direct approach, perhaps using 
the coordinates on the vortex moduli space introduced in [51].) As we have seen, the 
unique ground state of k vortices is given by (3.4). Here we are interested in excitations 
about this ground state which obey the linearisation of the constraint (3.2). These were 
described in [48]. They take the particularly simple form, 

5iZ = and 5i<p = 0 with l — l,...,k (B.l) 
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There are no higher excitations because {Z\) k = 0. We then describe the most general 
excitation as 

n 

Z{t) = Z, + Y J Ci(t)Zl l ~ l 
1=1 

with the complex coefficients q providing coordinates on the moduli space in the neigh¬ 
bourhood of the ground state. We work in a gauge with vanishing U(k) gauge field. 
Then the kinetic terms in (3.1) are 

Ski„ = j dt Tr \V t Z\ 2 + \V t <p\ 2 = J dt (Tr Z^Z^ 1 ) c*c n 

Meanwhile, the potential term from (3.1) is 

•S'pot = - j dt Tr | [a, Z] + uZ\ 2 + tp^a 2 tp = - j dt (Tr Z l ^~ 1 Zl n ~ 1 ) uj 2 l 2 c*c n 

where we’ve used (3.3) and the expression for o given below (3.4). We’re left having 
to compute the traces 

Tr Z[ 7 l Zl n ~ l = e t 6 ln with 0 ; = k(k - 1)... (k - 1 + 1) 

The low-energy dynamics of the vortices can then be written as 

k „ 

S = / dt (c* ci - u 2 l 2 c?a) 

i=i d 

This can be interpreted as the action for a real, massless scalar field, propagating around 
the edge of the vortex disc, truncated to the lowest k Fourier modes. The zero mode 
is also projected out. If the disc has radius R ~ Vk/ev , the scalar field has velocity 
v = ojR. As k increases, the radius of the disc scales as Vk, which means that the 
density of modes scales as 1 /y/k. This suggests that there is a continuum d = 1 + 1 
dimensional limit as k —> oc. 

The calculation above closely follows an analogous calculation for vortices in a first 
order system [52], The difference is that the edge modes for relativistic vortices prop¬ 
agate both ways around the circle while the non-relativistic vortices described in [52] 
have only chiral excitations. 

Here we have described the dynamics of £7(1) vortices in a trap. The dynamics of 
U(N) vortices seems somewhat harder. In the presence of masses , there are isolated 
ground states for the vortices, with the confined monopoles now acting as instantons 
interpolating between these vacua. It would be interesting to understand if there is a 
corresponding boundary interpretation for the resulting dynamics. 



C. Appendix: An Index Theorem for Non-Abelian Kinks 

In this appendix, we derive the result of the index theorem (3.8). 

The general path of the calculation follows previous index theorems for kinks in gauge 
theories. The index theorem for kinks in Abelian gauge theories was first derived in 
[40] and for non-Abelian gauge theories, in the absence of the adjoint scalar Z, in [29]. 
Both of these, in turn, closely follow the index theorem of E. Weinberg for magnetic 
monopoles [38, 39]. As we will see, the presence of the adjoint scalar Z makes the 
computation somewhat more involved. 

We start by rewriting the kink equations (3.6) in slightly more useful notation which 
allows us to avoid keeping track of the flavour indices i. We define a k x N matrix p 
whose columns are pi, and also a mass matrix m = diag(mi,..., mjv). Then we can 
write the kink equations in the more convenient form 

V 3 tp = (cr - m R )(p 

V 3 Z = [a,Z]+ujZ (C.l) 

2V = g 2 {[Z, z f ] + pp ] - r) 

where the subscript R denotes how the operator acts: m R p = (pm. We count the 
number of zero modes of these equations, which means that we count the number of 
solutions to the linearised kink equations. Given a background solution to (C.l) we 
write perturbations as hatted variables so, for example, p —» p + p. The linearised 
equations are then 

d 3 p = (a + ia 3 )p + (cr + ia 3 - m R )p 

d 3 Z — [<j + ia 3 , Z] + [a + ia 3 , Z] + uZ 

d 3 a = i[a 3 , a\ + i[a 3l a] + g' 2 ([Z, Z f ] + [Z, Z f ] + pp f + pp f ) 

To this we add a background gauge fixing condition, 


d 3 a 3 = i[a, cr] + i[a 3 , a 3 \ - ig 2 ([Z, Z f ] - [Z, Z f ] + pp f - pp f ) 


This naturally combines with the third linearised equation. We defined the complex¬ 
valued adjoint U(k) object £ = a + ia 3 , where a 3 is the spatial component of the U[k) 
gauge field and these can be written in the unified form 

%l=[l,f , ] + 29 2 ([Z,Zt]+w t ) 
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The number of solutions to these equations should be independent of g 2 ; to simplify 
the calculation we set 2 g 2 = 1. We then write the adjoint of these equations as 



(i*\ 



{ d,-£ a -<f> L -Z a \ 

A 


= 0 , 

A = 

-4>l d-^ R + m 


Uv 



{ ~Zl d + ti-u) 


Once again, subscripts denote the action of the operator: a means adjoint, with £ a X = 
[£, X]\ L means left multiplication, with <p\X = 0t X. Our goal is to find the number 
of solutions to these equations or, equivalently, the number of zero eigenvalues of A. 


We can also construct the adjoint operator to A, namely 


A f 


-OL -Z a \ 

-<j>L -d-^R + m 

V -Zl -d + Z a -tuJ 


Here there is a very important point: it is not true that At has no zero modes. This is 
in contrast to most index theorems for BPS solitons. It means that, for some range of 
parameters, the index that we compute undercounts the number of zero modes of the 
configuration; we give some examples in the main text. 


We define the regularized index 


X (M 2 ) 


Tr 


M 2 

At A + M 2 


M 2 

AAt + M 2 


where the trace is a functional one. The index that we want is 


X = lim X(M 2 ) 

M—>0 

This counts the number of zero modes. 

A standard set of manipulations [40, 29, 38, 39] allows us to express the index in 
terms of a total divergence of a current, and hence as a boundary term, 

where 


J(x 3 ,M 2 ) = Tv(x 3 \j\x 3 ) 


(C.2) 
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and 


3 = A- 


A f - 


At A + M 2 AAt + M 2 

Hence the problem reduces to computing this trace in this particular limit. Conve¬ 
niently, we only need the trace in the asymptotic regime where (£, Z, (ft) solve the 
vortex equations (3.2) and (3.3). In particular, we can use £ = £t = a. Asymptotically, 
the two squared operators are identical and take the simple form 

/©i \ 


AA f 




= A f A 


£ 3 =d=oo 


©2 

Z&L 


<p\ Z a 

©3 


(C.3) 


where 


©i = -d 2 + o 2 a + <p L <p[ + Z a Z\ 
©2 = -< 9 2 + {<jr - mf + ip[<p L 
@3 — —d 2 + (o a — Uj ) 2 + z\z a 


Note, however, that the presence of the adjoint scalar Z means that the matrix (C.3) is 
not block-diagonal. This makes the computation mathematically trickier than that of 
[29], but ultimately these terms will not contribute to the index. Instead, the possibility 
of fluctuations in Z are responsible for the richer structure of the index theorem 


The traces in (C.2) separate into two terms, 

-o R + m 


Tr j = -2 Tr o a (©i + M 2 ) 1 + 2 Tr 


o a - UJ 


©2 + M 2 ip{z a 
Z& L 0 3 + M 2 


-l 


To make progress, we need to use a number of properties of these matrices when eval¬ 
uated on the vortex solutions (3.5). We begin with the first term. One can show that 
[oai VLVl + ZaZl] = 0. This means that [<r 0 , ©i] = 0, and so we may simultaneously di- 
agonalise these two expressions. Denote the eigenvalues of o a by A a , and the eigenvalues 
of (p L (p[ + Z 0 Z\ by \' a . Then 

K, 


Tr o a (©! + M 2 ) =]T Z 


a -a 2 + M 2 + A 2 + A' 

A straightforward computation then shows that 

o a X = \ a X ^ a a A t = -A a A" t 
(<PL<pl + Z a Z\)X = X' a X =► {<p L <p[ + Z a Zl)X t = Alxt 


(C.4) 


As a result, we find that all terms in (C.4) pair up and annihilate each other. The first 
term in Tr j vanishes. 
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We can treat the second term in Tr j similarly. We observe that 


where 


@2 + M 2 (p\Z a \ 
Z&L ©3 + M 2 ) 


-d 2 + M 2 + Q 2 + R 


Q 


-a R + m 



R = 


PlVl P L Z a 
ZlfL zlz« 


satisfy [Q, R] = 0. Denote the simultaneous eigenvalues of Q, R as pa and p ' A . The 
number of solutions must be independent of simultaneously rescaling (</?, Z, y/r) and 
so cannot depend on the values of the p A . Technically this arises because, for each 
non-zero eigenvalue of R, the corresponding eigenspace has vanishing Q trace. This 
means that we can process by setting p! A = 0 to arrive at the following expression: 


J 



V 


Tr j 


d = £ 


Pa 


V M 2 + Pa 


s S n (^) 

A 


We’re left having to compute the signs of the eigenvalues of — cr R + m and <J a — u. 
Because a takes the diagonal form (3.5), the eigenvalues of — a R + m are simple: 

E ( m J _ m i + (q —!) w) 

i 

The eigenvalues of cr a — co are the differences of all pairs of diagonal elements of a 
(because cr a is in the adjoint representation). This gives 


hi kj 

££££ sgn (rrij — rrii + (q — p — 1) u>) 

i j q=\ p — 1 


Combining the two, we have our final expression for the index 


X = - 
2 


ki -1 kj 

£££ sign(mj — rrii + (q — p)u>) 

i,j q =0 p =0 




x^=—oo 


kj — 1 


££ sign(mj — rrii + (q — kj)u >) 
L i,j <2=0 


#—+00 


x' J =—oo 


where we’ve retained the first, more complicated version of the sum to illustrate the 
symmetry in the counting. By making use of the anti-symmetry in (i,j) we arrive at 
the second line. This is the result (3.8) in the main text. 
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